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Abstract 

A unique set of micromechanics equations Is presented for high temperature 
metal matrix composites. The set Includes expressions to predict mechanical 
properties, thermal properties, and constituent microstresses for the 
unidirectional fiber reinforced ply. The equations are derived based on a 
mechanics of materials formulation assuming a square array unit cell model of 
a single fiber, surrounding matrix and an Interphase to account for the 
chemical reaction which commonly occurs between fiber and matrix. A 
preliminary validation of the equations was performed using three-dimensional 
finite element analysis. The results demonstrate excellent agreement between 
properties predicted using the micromechanics equations and properties 
simulated by the finite element analyses. Implementation of the micromechanics 
equations as part of an Integrated computational capability for nonlinear 
structural analysis of high temperature multilayered fiber composites Is 
Illustrated. 

Key Words: Metal matrix composites; Composite micromechanics; Mechanical 
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Introduction 


The mechanical performance and structural Integrity of fiber reinforced 
metal matrix composites are ultimately governed by the behavior of the 
constituent materials at a mlcromechanlstlc level. In general, the Individual 
constituents behave quite differently relative to one another. Moreover, 
behavior of the constituents Is dynamic, particularly In high temperature 
applications, due to the various nonllnearltles associated with, for example: 
(1) large local stress excursions, (2) temperature-dependent material 
properties, (3) time-dependent effects, and (4) constituent chemical reaction. 

In the structural analysis of metal matrix composites, then. It Is 
Important to be able to describe and track this mlcromechanlstlc constituent 
behavior. Available methods for this purpose are limited. For example, 
techniques such as finite element analysis can, In principle, be applied 
directly with the constituents modeled discretely. It becomes obvious, 
however, that for complex structures the resources (manpower and computer) 
necessary to define, conduct and Interpret such an analysis are prohibitive. 
Another approach is to employ composite micromechanics theory and derive 
simplified relationships which describe the three-dimensional anisotropic 
behavior of the simple composite (e.g., unidirectional ply). The latter 
approach has been taken as part of a comprehensive research program to develop 
effective computational mechanics methodologies for high temperature 
multilayered fiber composite structures. 

As an essential part of the above-mentioned program, a unique set of 
micromechanics equations has been derived for high temperature metal matrix 
composites. The set comprises closed-form expressions to predict equivalent 
"pseudo homogeneous" properties for the unidirectional fiber reinforced ply. 
Including: (1) mechanical properties - moduli, Poisson's ratios, and uniaxial 
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strengths; (2) thermal properties - conductivities, coefficients of expansion, 
and heat capacity; and (3) constituent microstresses. 

The micromechanics equations presented here are derived based on a 
mechanics of materials formulation assuming a square array unit cell model of 
a single fiber, surrounding matrix and an Interphase to account for the 
chemical reaction which commonly occurs between fiber and matrix. The basis 
of the formulation Is summarized as part of the discussion below. 

Concurrent with the derivation of equations, a study was conducted using 
three-dimensional finite element analysis. The purpose of the study was to 
assess the validity of the mechanics of materials formulation, In general, and 
to Investigate the accuracy of the micromechanics equations for a specific 
composite material system. Results from this study are presented also as part 
of the discussion below. 

Finally, a demonstration of the utility of this unique set of 
micromechanics equations Is provided by Illustrating their use as part of an 
Integrated computational capability for the nonlinear structural analysis of 
high temperature multilayered fiber composites. A few typical results are 
presented from the stress analysis of a hypothetical tungsten fiber reinforced 
superalloy turbine airfoil. 

Composite Micromechanics Theory 

Composite micromechanics theory refers to the collection of physical 
principles, mathematical models, assumptions and approximations employed to 
relate the behavior of a simple composite unit (e.g., lamina or ply) to the 
behavior of Its Individual constituents. For example, a variety of approaches 
have been used In the past to predict equivalent thermoelastic material 
properties of unidirectional fiber composites [1-6]. More recently, simple 
equations have been derived [7,8] to predict mechanical, thermal, and strength 
properties for resin matrix composites using a mechanics of materials 
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formulation. A similar approach was taken to derive the set of micromechanics 
equations presented here for high temperature metal matrix composites. 

The formal procedure of composite micromechanics theory relies on the 
principles of solid mechanics, thermodynamics, etc., at different levels of 
mathematical sophistication, together with certain assumptions (consistent 
with the physical situation) and approximations. In the approach taken here, 
application Is made of the principles of displacement compatibility and force 
equilibrium as defined In elementary mechanlcs-of-materlals theory and 
Fourier's law for heat conduction from thermodynamics. In addition, the 
assumptions are made that: (1) fibers are continuous and parallel; (2) 
properties of all fibers are Identical; and (3) complete bonding exists between 
constituents. No restrictions need be placed on the constitutive behavior or 
Isotropy of the Individual constituent materials. For generality, constituent 
material behavior can be taken as thermoviscoplastic, anisotropic, and 
three-dimensional. It Is Implied by this that the Individual constituent 
material histories can be tracked Independently as a function of time and 
represented as an Instantaneous stress/straln state. 

The periodic structure of a unidirectional metal matrix composite (ply) 

Is approximated here by a square array unit cell model. The geometry of the 
model Is Illustrated In Fig. 1. It should be noted that the Interphase growth 
Is assumed to result from the degradation of fiber material and thus propagates 
Inward causing a continuous decrease of the current (Intact) fiber diameter (D) 
from the original (virgin) fiber diameter (D q ). With the existence of the 
Interphase, three subregions (A,B,C) are distinguished to characterize the 
Intralaminar (through-the-thlckness) nonuniformity of the constituent (matrix 
and Interphase) microstresses and material properties. 

The definition of ply properties Is with respect to the ply material 
coordinate system which Is depicted In Fig. 2. The common terminology 
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associated with each of the coordinate axis directions Is also Illustrated on 
the ply schematic. The micromechanics equations presented here are derived for 
the special case of a transversely Isotropic (Isotropic In the 2-3 plane) ply 
allowing for transversely Isotropic constituents. 

Composite Micromechanics Equations 

The micromechanics equations to predict ply equivalent mechanical 
properties are summarized In Fig. 3. Included are expressions for normal 
(extenslonal) moduli (E all , E l22 ) * shear moduli (G^, G l23 ), ancl Po1sson ' s 
ratios (v^.| 2 * v a23^’ In ex P ress ^ ons k represents constituent original 
volume fraction (values prior to any Interphase growth) and the subscripts f, 
m, d, and 8. denote fiber, matrix, Interphase, and ply quantity, respectively. 
The volume fraction of Interphase Is expressed In terms of the fiber original 
volume fraction and the virgin and Intact (In situ) fiber diameters. 

The equations for moduli are derived with modulus taken In the general 
context as simply the derivative of stress with respect to strain. As such, 
the expressions are applicable to the prediction of Instantaneous or tangent 
moduli as well as elastic moduli. It should be noted that the expressions for 
transverse moduli do not account for the longitudinal Poisson restraining 
effect that the fiber Imparts on the matrix. The restrained matrix effect Is 
considered here to be negligible for metal matrix composites. 

The effect Is generally more significant In resin matrix composites, for 
example, where the fiber/matrix relative stiffness ratio Is much greater. 

The ply equivalent thermal properties are predicted by the micromechanics 
equations summarized In Fig. 4. Included are expressions for heat capacity 
(C^), thermal conductivities (K^, K a22 ), an( * thermal expansion coefficients 
(a.,,, a i22 ). In the ex P ress1on for hefl t capacity the symbol p represents 
density. 
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The ply In-plane uniaxial strengths are predicted by the micromechanics 
equations summarized In Figs. 5 and 6. Included are expressions for tensile 
strength (S lllT , S^^), com P ress1ve strength ( s tl c . S t22C ). and 
Intralaminar shear strength Eac ^ oE p ^ y s ^ ren 9ths Is 

associated with a specific failure mode, as Illustrated by the schematics In 
Fig. 7. In the case of longitudinal compressive strength, four different 
failure modes are considered. The four expressions In Fig. 5 for 
correspond, respectively, to the four failure modes as follows; fiber 
compression mode, matrix compression mode, delamlnatlon/spllttlng mode, and 
fiber mlcrobuckllng mode. A more comprehensive treatment of micromechanics 
strength theories Is given by Chamls [9]. 

The expressions to predict the thermomechanical microstress distribution 

In the ply constituents are summarized In Figs. 8 to 10. Included are 

expressions for fiber microstresses (o^, a^,,, a fl2’ °f23^ ^erphase 

6 C 8 C B C 

microstresses (o^, °d22* 0 dl2’ °d2p and ma * r1x mlcrostresses 

°m22 ,C ’ °ml2 ,C * °m23 C ^' * n ex P ress ^ ons AT represents an Incremental 
change In temperature and the superscripts A, B, and C denote the 
Intralaminar subregions Illustrated In the accompanying schematics. It should 
be noted that these expressions for constituent microstresses are based on 
uniaxial behavior, l.e., they do not Incorporate any Poisson contributions. 

The systematic procedure for deriving the micromechanics equations 
summarized above Is explicitly demonstrated In the Appendix with the 
derivations for normal moduli (E^ and E^). Der1vat1ons of the other 
equations are omitted here solely for the sake of brevity. The selection of E^ 
and for demonstration purposes was based on the authors' judgment that 

their derivations are sufficiently representative to adequately demonstrate the 
formal procedure. 
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Mlcromechanlcs/Flnlte Element Validation 

In order to Investigate the validity of the mechanics of materials 
formulation and assess the accuracy of the equations derived therefrom, a 
preliminary study was conducted using three-dimensional finite element 
analysis. The objective of the study was to compare the equivalent ply 

properties (E lir E l2 2’ G ju 2’ G i23’ Yl2’ V 123’ “ill’ a t22 ) P red1cted 
by the micromechanics equations with the average "pseudo homogeneous" ply 

properties simulated In the finite element analyses. 

To conduct the analyses, a discrete model of the square array unit cell 
was constructed, as shown In Fig. 11, from Isoparametric solid finite elements. 
The composite material system assumed for this study Involved a thorlated 
tungsten (W-1.5Th0 2 ) fiber embedded In an Iron-base superalloy (Fe-25Cr-4Al-l Y) 
matrix. Properties for the Interphase were taken to be a simple average of the 
fiber and matrix properties. 

The analyses entailed simulations of Idealized modes of deformation such 
as simple elongation, pure shear, and unconstrained thermal expansion. These 
were achieved through the judicious application of the loading/boundary 
conditions on the model. The appropriate simple expressions from elementary 
mechanics of materials theory (see Fig. 11) were then applied In conjunction 
with the nodal displacement/force results of the finite element analyses to 
compute the simulated average properties of the discrete model as a "pseudo 
homogeneous" unit. 

Results of the study are summarized In Table 1 which gives the ratios of 
property values determined from the micromechanics equations (P uirn ) and by 
finite element simulation (Pp EH ). As can be seen, excellent agreement was 
achieved overall. These results Indicate that the mechanics of materials 
formulation Is an effective approach to the micromechanical modeling of metal 
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matrix composites. It Is recognized, however, that additional Investigation, 
both analytical and experimental, would be prudent before any final conclusions 
are made regarding the specific accuracy of these micromechanics equations. 

Application of Micromechanics Equations 
The primary Impetus In deriving the set of micromechanics equations 
presented here was for Implementation as part of an Integrated computational 
capability for the nonlinear analysis of high temperature multilayered fiber 
composites [10]. This particular utilization of the equations Is demonstrated 
here with a few typical results taken from the nonlinear (quasi-static) stress 
analysis of a hypothetical turbine blade (airfoil only) model. The 
Incremental/lteratlve analysis was conducted to Investigate the thermally 
Induced residual stresses developed during the cool-down transient of a typical 
fabrication process. 

The airfoil Is a hollow thin shell structure of constant thickness with 
walls comprising a four-ply [±45] $ laminate based on W-l.SThO^ fiber 
reinforced Fe-25Cr-4Al-l Y at a fiber volume fraction of 0.50. Since the 
purpose here Is merely to Illustrate the types of Information provided by the 
micromechanics equations In this particular Implementation, further details of 
the airfoil model and analysis are omitted. 

Two examples of ply mechanical property predictions are given In Figs. 12 
and 13 which show the variation during the cool-down transient of constituent 
and ply longitudinal and transverse moduli, respectively. The ply moduli are 
computed from the corresponding micromechanics equations. The results In 
Fig. 12 reflect the rule-of-mlxtures relationship expressed by the equation for 
while the results In Fig. 13 Illustrate the dominance of the matrix 
modulus on the value for £^ 2 ' 
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The development of residual stresses during the cool-down transient Is 
Illustrated In Figs. 14 and 15. The results are for the longitudinal and 
transverse normal components, respectively, of ply stress and constituent 
microstresses. The microstresses are computed from the corresponding 
micromechanics equations. The points to be noted from these results are the 
relative magnitudes and sense (tensile or compressive) of the constituent 
microstresses. In Fig. 14, for example, the opposite sense of the fiber and 
matrix microstresses results from the difference In thermal expansion 
coefficients between the two materials. The results In Fig. 15 Illustrate the 
significant through-the-thlckness nonuniformity of the matrix and Interphase 
microstresses, as characterized In the different Intralaminar subregions 
(A,B,C) . 

From just the few examples given, the utility of the micromechanics 
equations becomes more apparent. Considering the results of microstress 
distribution, for example, It becomes Intuitively more clear how material 
failures might occur at a local level and prompt the Initiation of a flaw. 

This type of Information provides an Insight Into the behavior of composites 
at a mlcromechanlstlc level which undoubtedly Influences their performance and 
Integrity In a structural application. 

Summary 

The set of micromechanics equations presented here for high temperature 
metal matrix composites Includes expressions to predict the mechanical 
properties, thermal properties, and constituent microstress distribution for a 
unidirectional fiber reinforced ply. The equations Incorporate an Interphase 
region at the fiber/matrix boundary In order to account for the chemical 
reaction which commonly occurs In high temperature applications of these 
composites. The basis of the mechanics of materials formulation from which the 
equations are derived Is described. The formulation Is shown to be a valid 
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and effective approach to micromechanical modeling of metal matrix composites, 
supported by the favorable results achieved In a comparison with 
three-dimensional finite element analysis. The utility of the micromechanics 
equations as part of an Integrated composite structural analysis capability Is 
Illustrated with examples taken from the nonlinear stress analysis of a turbine 
airfoil. The results demonstrate the ability to describe and track behavior 
at a mlcromechanlstlc level which Impacts the performance and Integrity of 
these composites In structural applications. 
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Appendix 

In order to demonstrate the formal procedure Involved In the application 
of composite micromechanics theory, derivations of the equations for ply normal 
moduli (E all and E fc22 ) are explicitly developed below. The particular approach 
taken here relies on the principles of force equilibrium and displacement 
compatibility as defined from elementary mechanlcs-of-materlals theory. 

Longitudinal Normal Modulus 

Consider the square array unit cell model (see Fig. 1) subjected to a 
uniaxial load In the longitudinal direction (see Fig. 2). The equivalent 
composite (ply) load Is defined from force equilibrium to be the sum of the 
constituent loads as follows: 


P l = P f tP d + P m H) 

In the Integrated average sense, Eq. (1) Is rewritten as 

o. A. = O.A.. + o. A. + o A (2) 

aaffddmm . ' ' 

where A represents cross-sectional area. Dividing through by A^ and 

noting that because of a common longitudinal dimension the resulting area 

ratios are equivalent to actual volume fractions, Eq. (2) reduces to 

°1 ■ °f k f * °d k d * °m k l < 3 > 

Because compatibility of longitudinal displacement requires equal strains for 
the composite and constituents (c^ = = e m ), Eq. (3) can be 

differentiated with respect to strain to give 


'do. 


do 


£ K‘ 


'do , 


do_ 


k d + 


m 


dc 


m 


(4) 


The quantities (do/de) represent the slopes of the corresponding stress- 
strain curves for the composite and constituents and In this context define 
Instantaneous or "tangent" moduli. Hence, Eq. (4) becomes 


E. = E,k, + E .k . + E k 
9. f f d d mm 


(5) 
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Expressing actual volume fractions In terms of original fiber and matrix volume 
fractions (before Interphase growth) and original and Intact fiber diameters, 
Eq. (5) Is rewritten as 



kE 
m m 


( 6 ) 


Equation (6) Is the desired form and Is the same as that given In Fig. 3. 

Transverse Normal Modulus 

Consider the square array unit cell model again except that the fiber and 
Interphase are of equivalent square cross-section such that linear dimensions 
(In the plane of cross-section) can be defined as follows: 


a f 




(7) 


and 

Sf . a f , s d . a d - a f . s, - a, - a d> a, ■ a, (8) 

Assume a uniaxial load In the transverse direction and neglect Poisson effects. 
For subregion C displacement compatibility yields 

Vs. = s f c f + s d c d + s m c m (9) 

and force equilibrium results In equal stresses for the composite and 

constituents (a. = a, = a . = a ) . Hence, eq. (9) can be differentiated 
8, f d m 

with respect to stress to give 



The quantities (de/do) represent reciprocals of the slopes of the 
corresponding stress-strain curves for the composite and constituents and In 
the same context as before define reciprocals of Instantaneous or "tangent" 
moduli. Hence, with some rearranging Eq. (10) becomes 
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m 

A\ „ N\ (la\ t fr\ /EbY' 

AV \ s l/ \ E d/ \ S lA E f/. 


Substituting the definitions In Eqs. (7) and (8) Into Eq. (11) and rearranging 
gives 


which defines an equivalent modulus for subregion C. The equivalent modulus 
for subregion B Is deduced from Eq. (12) by letting 0/0 Q equal unity. The 
result Is 


1 -V*7 1 - 


The equivalent modulus for subregion A Is simply the matrix modulus or 

£ f - E m <1,; 

The ply transverse modulus (^^22^ * ^ en » defined by assuming that 
subregions A, B, and C act as parallel elements when subjected to a transverse 
load. This Is analogous to the case for E^ where the constituents are 
assumed to act In parallel. Hence, from Eq. (5) It Is deduced that 

E l s t - E l s f * E l s d * E 4 s m < 15 > 

Dividing through by sj,, substituting the definitions from Eqs. (7) and (8) 

and the results from Eqs. (13) through (15), and rearranging gives 


VM 1 - d7 


V =E m (1 - k f> + 


f ) + — r / nr + r — i — vn — tx /m 

1 -v^f 1 - rl 1 -V*f 1 - 1 - r r - rlc 

L \ f/J L \ o/ \ d/ \ o/ \ f/J 


Equation (16) Is the desired form and Is the same as that given In Fig. 3. 
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TABLE 1. - MICROMECHANICS/ 
FINITE ELEMENT VALIDATION; 
COMPARISON OF PROPERTY 


PREDICTIONS/ 

SIMULATIONS 


Property 

P /P 
MEQ / FEM 

E 8.11 

1.00 

e p22 

1.01 

G a.i 2 

.96 

g 123 

.98 

m 2 

1.00 

v &23 

1 .08 

“111 

.99 

a a22 

1 .15 


P MEQ _ Property 
predicted by 
micromechanics 
equation. 

PpEM - Property 
simulated by finite 
element analysis. 
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Figure 1. - Micromechanics model; square array unit cell. 
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Figure 2. - Unidirectional composite (ply) material coordinate 
system. 



THROUGH-THE-THICKNESS 

! t 



LONGITUDINAL 




Figure 3. - Micromechanics equations; ply mechanical properties. 



H 

K m 

K iZ2' 

“ 111 ' 
a t22 ’ 


THROUGH-THE-THICKNESS 





Figure 4 -Micromechanics equations; ply thermal properties. 
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Fiqure 5. -Micromechanics equations; ply uniaxial strengths, longitudinal. 
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Figure 6. -Micromechanics equations; ply uniaxial strengths, transverse and shear. 



(a) Longitudinal 
tension. 



(b) Longitudinal compression. 


r.bi.b.oidJo.'O-. 

:OlVo* 6 :d; 0 . a 
^ o. b- o?'d jrt oj- 

.fi-Cj /J-.fT. 

,o *o 1 q vo-f o 'ib. 

;'o io -o jo- b u). ? 

. 0: rv?o:d 


iO- (V h i>. n 




I > 




6 o,u6.0 o.o yjt. 


b"o a 'p’.o- b J 
T 'o~li r C>' T <>r rt-b*^q"b *° •bjlbi'f 
p'VFtoVq. 6 '6: q^o ;p ‘VjA Ai. 

.VAW&5?« 


II 


(c) Transverse tension. (d) Transverse compression. (e) Intralaminar shear. 

Figure 7. - In-plane failure modes for unidirectional ply. 
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Figure 8. - Micromechanics equations; fiber microstresses. 
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Figure 9. - Micromechanics equations; interphase microstresses. 
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Figure 10. - Micromechanics equations; matrix microstresses. 




Figure 11. -Micromechanics/finite element validation; finite 
element model and simple mechanics of materials expressions 
for idealized deformation modes. 
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Figure 12. - Fabrication cool-down transient; 
variation of for constituents and ply. 
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Figure 13. - Fabrication cool-down transient; 
variation of E 22 for constituents and ply. 
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Figure 14. - Fabrication cool-down transient; 
induced residual stress (o^) for ply and 
constituents. 
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Figure 15.- Fabrication cool-down transient; 
induced residual stress (oW for ply and 
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The following corrections apply to the appendix and occur on page 13: 

1. The denominator of Equation (12) should read as follows: 



2. The first full sentence after Equation (12) should read as follows: 

"The equivalent modulus for subregion B Is deduced from Eq. (12) by 
letting D/D 0 equal zero." 


3. The denominator of Equation (13) should read as follows: 



4. The sentence after Equation (15) should read as follows: 

"Dividing through by Sj,, substituting the definitions from 
Eqs. (7) and (8) and the results from Eqs. (12) through (14), and 
rearranging gives" 

5. The denominator of the second term Inside the braces on the right-hand 
side of Equation (16) should read as follows: 
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